Rapid laser-free ion cooling by controlled collision 
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I propose a method to transfer the axial motional excitation of a hot ion to a coolant ion with 
possibly different mass by precisely controlling the ion separation and the local trapping potentials 
during ion collision. The whole cooling process can be conducted diabatically, involving only a few 
oscillation periods of the harmonic trap. With sufficient coolant ions pre-prepared, this method can 
rapidly re-cool ion qubits in quantum information processing without applying lengthy laser cooling. 



Ion trap system is an auspicious implementation of 
quantum computers due to its excellent controllability 
and stability 1 , J,] . Various building blocks of a quantum 
computer [3|, Ij] , such as fast and precise quantum gates 
0j J long time storage of quantum information [tI , and 
high fidelity readout [8] , have been demonstrated exper- 
imentally. In some scalable ion trap quantum computer 
architecture, ion qubits are shuttled across traps for dif- 
ferent operations iS]. During computation, the qubits 
would be heat up by various mechanisms, such as fluc- 
tuations of trapping potential, imprecise transportation 
of ion qubits, and momentum gain in fluorescence read- 
out, which will diminish the performance of subsequent 
logic operations. Ion qubits are usually re-cooled by sym- 
pathetic cooling Typical duration of current state 
sympathetic cooling is at millisecond range which 



is at least two order of magnitude longer than the time 
scale of other quantum operations (13] . The time bot- 
tleneck of sympathetic cooling should be resolved for a 
faster ion trap quantum computer for better preservation 
of quantum coherence and higher computational power. 

Here I present an alternative approach that takes only 
a few trapping oscillation periods to cool an ion. The 
idea is to interchange the axial motional state of an ion 
with that of a coolant ion through a phonon beam splitter 
[l3| . which can be implemented by precisely controlling 
the ion separation and the local potentials. A ground mo- 
tional state can then be transferred from the coolant ion 
to the hot ion through their mutual Coulomb interaction, 
without the aid of laser. I will outline the procedure to 
devise the time variation of the trapping parameters for 
the cooling process. Although similar idea has been pro- 
posed and experimentally demonstrated in 14 , , their 
operations are adiabatic in order to avoid parametric ex- 
citations. The speed of my method can be beyond the 
adiabatic limit as the problem of parametric excitations 
can be solved by using the dynamic invariant formalism 
of time dependent harmonic oscillators [li^. My cooling 
operation can be boosted up by bringing two ions closer 
for stronger Coulomb interaction. 



I consider two ions are radially tightly trapped but 
axially weakly trapped by a double well potential, which 
can be implemented by microfabricated surface traps 17|, 
18|. The Hamiltonian is given by 
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where ion 1 is the ion to be cooled and ion 2 is the coolant. 
The ions can be different in species, i.e. mi ^ TO2. The 
trapping parameters lui, uj2, Ri, and R2, are related to 
the local potential experienced by the ions and are as- 
sumed to be tuneable. 

By separating the classical and the quantum attributes 
of the ions' motion [l3|, i.e. Xi = Xi(t) + qi for Pi = 
Piit) + T^i for i = 1,2, the classical equation of motion 
can be extracted as 



rriiXiit) = -mtujf{t) {x.i{t) - Ri{t)) 
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where the ion separation r(t) = xi{t) — X2(t) > 0. 
The quantum attributes of motion follow the Schrodinger 
equation 
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The last round bracket denotes a Taylor expansion of the 
Coulomb potential, which is valid if the mean quantum 
attributes of position, {qi), are much smaller than the 
minimum ion separation. If only the zeroth order of the 
expansion series is considered, the Hamiltonian is of sec- 
ond order in quadrature operators. Then the effect of 
the evolution operator, /7(t), is equivalent to applying a 
gen eral two-mode squeezing operator on the initial state 

The aim of this paper is to determine the initial state of 
ion 2 and the time variation of cj's and i?'s such that ion 
1 will result in the motional ground state. Eq. ([3]) can 
be easily solved in a time independent squeezed basis. 
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i.e. IV') = S\^), where S' is a squeezing operator acting 
on ion 2 only and 7r2 — S^Tr2S = \J nil j mVk'^^ (72 — > 
S'^<l2S — \J mX I ■m2c(2 (the operators in the squeezed ba- 
sis are primed.) I introduce the centre of mass mode 
(-1- mode), a'^ — (a'j^ + al^)!^/!., and the breathing mode 
(- mode), d'_ = (o!-^ — a!<2)l\/l, where a is the annihi- 
lation operator defined by the quadrature operators in 
respective basis, viz. tt — —i^/m\hjjQ/2{a — a)) and 
q = \Jh/2miijjQ{a + a^)\ ojq is an arbitrary constant pa- 
rameter that can be taken as the frequency of the storage 
trap of both ions. The state |$) follows the equation 



which satisfies 



h-{t) + Ujl{t)b-{t) 
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Let the collision runs from t = —T/2 to T/2. For our 
purpose of a phonon beam splitter, the auxiliary func- 
tion should satisfy 6_(t < -T/2) 1 and 6_(t > 
T/2) ^ 1 for minimal parametric excitation, and 6- = 
I-T^/2 ^o/b'^{t)d.t = TT + ljqT for complete motional state 
transfer. Constructing a a;„(<) for a 6_(i) that satis- 
fies all the boundary conditions is possible but not easy. 
Instead, w-(t) can be inferred from a conjectured 6-(t) 
^ that already satisfies the boundary conditions. 

With 0Ji{t) and u}2{t) are set to satisfy respective con- 



,(4) 



up to the second order of g's. The mode frequencies are 



4(t) = 2KW 



c|(0) 



(5) 

Eq. Q is a two-mode coupled time dependent har- 
monic oscillator, of which the general analytic solution 
is not known yet. However, the system can be decou- 
pled into two single mode harmonic oscillators, i.e. the 
coefficient of q'_^q'_ vanishes, if wc know the classical tra- 
jectory of the ions and set the local potential of ion 2 as 
Lo^it) = ujlit) + (1/mi - l/TO2)eV27reor3(t). 

The main idea of the cooling operation in this paper is 
as follow: ion 2 is initially prepared as the ground state 
in the squeezed basis, in other words a squeezed ground 
state, 5*10), in laboratory; then a phonon beam splitter 
in the squeezed basis is applied to transfer a ground state 
to ion 1 i.e. an operation that transforms — a[ 
up to some unimportant phases. 

A squeezed ground state can be constructed by, for 
example, varying the trapping potential 13!, l21| or ap- 
plying a parametric drive 2l|. Both methods have been 
demonstrated in experiments 22, 23l |. 

A phonon beam splitter can be implemented by con- 
trolling Lu{tya and i?(t)'s, such that the ions collide 
in a way that their mutual Coulomb interaction trans- 
form the annihilation operators as a'_^_ — ^ a'^e~*^+ and 
d'_ — > a'_e~^^~ [ll]. The phase difference should be set 
as 0_ — 0+ = TT for a complete state transfer. For sim- 
plicity, I set iL'+{t) — luq, and so a'^ — > d'.^e"*'^"* after 
some operation time t. This requires uJiit) to be tuned 
to satisfy LL>l{t) = - (1/mi - ^l/mim2)e^/27reor^(i). 

The time variation of can only squeeze and 

phase-shift the - mode. According to the dynamic in- 
variant formalism of time dependent harmonic oscilla- 
tors [l6^, the squeezing parameter and the shifted phase 
are determined by a real scalar auxiliary function, &-(i). 



straints, oj-{t) depends only on r{t) according to Eq. ([5|). 
Therefore the conjectured fo-(t) also specifies the classi- 
cal motion of the ions. The desired r{t) can be attained 
by tuning Ri{t) and i?2(i) according to Eq. Pro- 
ceeding this derivation requires an additional constraint, 
which could be xi{t) = —X2{t) or else for convenience. 

The boundary conditions that — > 1 before and 
after the collision specify the large initial and final ion 
separation. The initial and final speed of the ions are, 
however, difficult to be included as boundary conditions 
on the auxiliary function. The speeds can be obtained 
only by numerically integrating Eq. (jG]) and ([2|). 

Before the controlled collision, the ions are transported 
from storage traps or interaction zones by moving har- 
monic wells 9] . Likewise after the collision, ion 1 is trans- 
ported away to be stored or to proceed the next opera- 
tion, while ion 2 is transported away to be discarded or 
re-cooled. In these transportation stages, the Hamilto- 
nian of the moving harmonic wells is the same as Eq. ([I} . 
In the forward trip, the moving harmonic wells should ac- 
celerate the ions' classical speed to match that specified 
by b^{—T/2) at the position Xi{—T/2); in the backward 
trip, the wells decelerate the speed from that specified 
by 6_ {T/2) to whatever required in the next operation. 
For a smooth transition between the transportation and 
the collision stage, the trapping parameters i?'s and w's 
ought to be continuous at t — —T/2 and T/2. Because 
r/2) and 6_(T/2) have to be close to 1, wi and 
UJ2 can be tuned as ujq during the transition. Maintain- 
ing the oscillation frequency of the wells as ujq through- 
out the transportation stage may be experimentally con- 
venient but not necessary, as the parametric excitation 
in this stage is local and can be undone by subsequent 
squeezing operations. The continuity of i?'s, ions' posi- 
tion and, speed can be incorporated as constraints when 
constructing the trajectories of the moving potential well. 
The ions can be transported diabatically with arbitrary 
speed if the harmonic wells are accurate enough. Such 
trajectories can be obtained by inverse-engineering or 
bang-bang method [241, or further optimised for addi- 
tional constraints j25„ ,26.] . 
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I now discuss about possible errors and the speed limit 
of the the cooling process. Because the ions are moved 
diabatically in the transportation stage, the speed of this 
process is not limited at the range of the trapping oscil- 
lation period [25|. Distortion to the quantum states can 
be caused by the anharmonicity of the trapping poten- 
tial and the imprecise control of the trapping parameters. 
As demonstrated in recent experiments, ions in MHz trap 
can be transported in a few fis without significant heat- 
Errors in the transportation stage are not 
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likely to reduce the performance of the cooling process 
at current level of technology. 

In the collision stage, errors may come from the infi- 
delity of ion 2's initial state, the inaccurate classical po- 
sition and speed of the ions at the stage transition, the 
anharmonicity of the Coulomb and the trapping poten- 
tial, and the fluctuation of the trapping potential. 

If ion 2 is prepared in a state pin instead of S'|0), then 
ion 1 will acquire the state pi^S, and thus the final ex- 
citation of ion 1 will be TT{a''aS'' pi^S}. I note that if the 
initial state is an evolving squeezed state, i.e. U2{t)S\0) 
where C/2(t) is the evolution operator of ion 2 in a har- 
monic oscillator with UJ2 = wq, it is not difficult to see 
that ion 1 still acquires a ground state after the process 
with an unimportant evolving phase. 

The inaccurate classical position and speed of the ions 
can be treated as an unknown quantum displacement, 
which effectively increases the excitation of the initial 
ions. In the forward trip, the displacement on ion 1 does 
not cause final excitation because it will be transferred to 
ion 2, while the displacement on ion 2 will be transferred 
to, and thus heat up, ion 1. Likewise in the backward 
trip, the displacement on ion 1 will cause final excitation 
while ion 2's state is not important. Nevertheless, the 
heating can be removed by subsequent displacement and 
squeezing operators if the displacement is known. For 
better performance of the cooling, the classical position 
and speed of the ions should be cognisant of at quantum 
level precision. 

To perturbatively study the performance of the cool- 
ing process under anharmonicity and fluctuations of po- 
tentials, I numerically simulate the collision between a 
""^Ca^ ion and a ^''Mg^ coolant ion. The simulated col- 
lisions are specified by 6_(t) = l/(e-*'/V%/A?7ri + 1) 
with fc = 2 (case I), fc = 3 (case II), and b-{t) = 
l/v/7exp(-t4/A:') + 1 (case III) with fc' = 8 and 7 sa 
3.048 that satisfy the tt phase difference condition. The 
simulations run for T sa 9/wo, 10.6/ajOi S.S/wq for case 
I, II, III respectively, where the value of T is defined 
b y r(-r/2) = r(T/2) = 20^0; wq = 27r MHz; Ig = 
•y/ e^/47^eo"^ll^o ~ 4.45/im. The time variations of the 
auxiliary function, the ion separation, and the trapping 
parameters are shown in Fig. [TJ 

The anharmonicity of the Coulomb potential are the 
higher order terms in the Taylor expansion in Eq. ([3]). 
The anharmonic terms become significant if the mean 




FIG. 1: Time variations of (a) (b) r, (c) Ri (solid line) 
and R2 (dashed line) i27|, (d) uji (solid line) and uj2 (dashed 
line), for the b-{t) of case I (blue), II (red). III (brown). We 
have required xi{t) — —X2{t). 
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FIG. 2: Final motional excitation of ion 1 caused by the 
Coulomb anharmonicity when the two ions are initially pre- 
pared in |nin)iS'|0)2 and then undergo the controlled collision 
specified by case I (blue), II (red). III (brown). 



initial excitation (fii)in is higher or the ion separation 
is shorter. Although the anharmonic heating is seem- 
ingly more serious in a faster cooling process that the 
ions should be brought closer for a stronger Coulomb in- 
teraction, my numerical results show that the heating 
effect is more sensitive to the form of &-(i) rather than 
only the speed. In Fig. [2l the process of case I yields 
generally higher excitation than that of case HI although 
the former one is slower. This result could be under- 
stood from Fig. [IJb) that the minimum ion separation 
in case I is shorter than that of case HI. In either case, 
the final phonon number of ion 1 is below 10^'^ even it 
has more than 40 phonons initially. We believe a process 
with faster speed but lower Coulomb anharmonic heating 
can be realised by employing a further optimised b-{t). 

The double harmonic well I consider in Eq. ([T|) is 
an approximation to the applied potential. In prac- 



tice, the applied potential may contribute anharmonic- 
ity to the quantum attributes. The magnitude of 
the anharmonic terms highly depends on the config- 
uration of the experiment. For example, if the ap- 
plied potential is a fourth order polynomial of posi- 
tion, i.e. V{x) = J2'j=i^j^''' where Vj are real pa- 
rameters related to w's and i?'s, then the leading or- 
der anharmonic terms in Eq. ([3]) are third order of g's, 
viz. [{miiol — m2w|)/4r + (micjji?i — m2a;|i?2)/?'^]<?i + 
[{miLul — m2Cj|)/4r — {miUjlRi — m2a;fi?2)/''^]92- I have 
tried that these terms cause (ni) to be only a few times 
higher than that caused by the Coulomb anharmonicity. 
However, these third order terms can be suppressed if 
the applied potential is a higher order polynomial of x. 
Thus, the actual magnitude of the anharmonic heating 
effect is very sensitive to the experimental configuration. 

Fluctuation of the trapping parameters can cause un- 
certainties on the ions' classical speed and position, and 
parametric excitation on the motional states. As men- 
tioned before, an unexpected classical displacement can 
be treated as a quantum displacement error a, so the an- 
nihilation operator transforms as Si — > ai + a. The final 
excitation of ion 1 becomes Ti{a\aip} — >■ Tr{a\dip} + 
jap, where the terms with Trjoipi} can be neglected 
if we consider the initial state of ion 1 is symmetric 
among the quadratures, e.g. thermal state, as the state 
remains symmetric if it follows Eq. I numerically 

estimate the order of magnitude of ion I's final exci- 
tation if w's and i?'s experience Gaussian fiuctuations 
with different frequencies and magnitudes. The results 
for case I are plotted in Fig. |3l As expected, heat- 
ing is more serious when the frequency of fluctuation is 
comparable to ujq, while high frequency noise heats in- 
significantly. A 10^'* accuracy of w's is enough for final 
excitation to below 0.1 phonon in the frequency range 
of fluctuation considered. For the fluctuation of i?'s, an 
accuracy better than quantum level of displacement, i.e. 
{a = l\q\a = 1) = y^2h/miU!Q where |a = 1) is a co- 
herent state of a harmonic oscillator with mass mi and 
frequency uq, should be needed for less than 1 phonon of 
final excitation. 

For the parametric excitation, the fiuctuation breaks 
the decoupling between -I- and - mode, but the Hamilto- 
nian in Eq. ^ remains quadratic in ^-'s and g's that guar- 
antees the annihilation operators transform linearly [l^ . 
So we can write U^t)aiU{t) = Aai + Bal + C 80,28^ + 
DSdlS'f for some complex A,B,C,D If the initial 
state of ion 2 is 5*10), and that of ion 1 is symmetric 
among the quadratures and has a mean phonon num- 
ber (rii)in, then the final motional excitation becomes 
Tr{d\aip} = + \B\^){ni)in + + \D\^. The sim- 
ulated results in Fig. [3] show that both \A\'^ + and 
are small parameters that the heating by para- 
metric excitation is comparable to that by displacement 
error if (fii)in is at the order of 10^. We conclude that 




FIG. 3: (a) and (b) show the heating of displacement un- 
certainties on ion 1, jap. (c) and (d) show the parametric 
excitation parameters \A\'^ + \B\^ (solid line) and |_Bp + \Df 
(dashed line). In (a) and (c), the w's are multiplied by a 
fluctuation, i.e. uji LJi{l + Suji), where Su's obey Gaussian 
distribution with standard deviation a = 10~^ (purple) and 
a = 10"^ (green). In (b) and (d), the _R's are added by a 
fluctuation, i.e. Ri ^ Ri + SRi, where SR's obey Gaussian 
distribution with width 2^2h/miLJo (purple) and ^J^KfnvJjJi^ 
(green). The Gaussian perturbation in both cases is time 
uncorrelated and updated for every time r in the numerical 
simulation. All the plotted data are averaged over 100 runs. 



the problem of parametric excitation is less important 
than the uncertainties of classical displacement in ion 
trap quantum information experiments. 

I now discuss the application of the cooling scheme. It 
is useful in rapidly re-cooling ion qubits during quantum 
computation. With sufficient coolant ions pre-cooled, a 
new one can be loaded in and collide with the qubit ion in 
each round of cooling [2^. Thus, the operational speed 
of the ion trap quantum computer is not limited by the 
duration of laser cooling processes. The scheme can also 
rapidly cool a pair of ions with the same mass. The ions 
are first split into two separate traps. Controlled collision 
is then conducted on each trap to bring the ions to the 
ground state. Finally they are re-combined by reversing 
the heating-less diabatic ion separation process [l3| . 

In conclusion, I propose that the axial motional excita- 
tion of an ion can be rapidly removed by controlled col- 
liding with a coolant ion that may have a different mass. 
I have outlined a procedure to devise the time variation 
of the trapping parameters for the process. Because the 
ions are not required to interact adiabatically, the cooling 
process can be conducted within a few oscillation periods 
of the trapping potential. 

I thank Daniel James for useful discussions. I would 
like to acknowledge support from the NSERC CREATE 
Training Program in Nanoscience and Nanotechnology. 
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